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Abstract 

The joint probability distribution function (PDF) of the height and its gradients is derived for a 
zero tension d + 1-dimensional Kardar-Parisi-Zhang (KPZ) equation. It is proved that the height's 
PDF of zero tension KPZ equation shows lack of positivity after a finite time t c . The properties of 
zero tension KPZ equation and its differences with the case that it possess an infinitesimal surface 
tension is discussed. Also potential relation between the time scale t c and the singularity time scale 
t c ,w-fO of the KPZ equation with an infinitesimal surface tension is investigated. 
PACS: 05.45.-a, 68.35.Ja, 02.40.Xx. 
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1 Introduction 



Studying the morphology, formation and growth of interfaces has been one of the recent interesting fields 
of study because of its high technical and rich theoretical advantages. On account of the disorder nature 
embedded in the surface growth, stochastic differential equations have been used as a suitable tool for 
understanding the behavior of various growth processes. Such equations typically describe the interfaces 
at large length scales, which means that the short length scale details has been neglected in order to derive 
a continuum equation by focusing on the coarse grained properties. A great deal of recent theoretical 
modeling has been started with the work of Edward and Wilkinson [1] describing the dynamics of height 
fluctuations by a simple linear stochastic equation. By adding a new term proportional to the square 
of the height gradient, Kardar, Parisi and Zhang made an appropriate description for lateral interface 
growth [2]. The d + 1-dimensional forced KPZ equation is written as 

d ±-^Vh? = vV*h + f (1) 

where h(x, t) specifics the surface height at point x ( d-dimensional vector ) and a > 0. The force / is a 
zero mean, statistically homogeneous, white in time, Gaussian process which it's covariance would be 

</(x, t)/(x', t')) = 2D D(x - x')5(t - t') (2) 

Typically the spatial correlation of the forcing is considered to be a delta function, mimicking the short 
length correlation. Here the spatial correlation is considered as 

^-^ = ^^oxp(-^-i^) (3) 

where a is much less than the system size L, i.e. a « L, which represents a short range character for 
the random forcing. It is useful to rescale the KPZ equation as h! = h/h , r' = r/r and t' = t/t . If we 
let h = (-77 1 ) 1 / 2 an d to = 17 1 where r is a characteristic length, all of the parameters can be eliminated, 
except the coupling constant g = a J? . The limit g — > 00 (or zero tension limit, v — > 0), is known 
as the strong coupling limit [13,14]. Although originally, this equation appeared as a model for surface 
growth [3], it is mostly used today in polymer physics [4], Burgers turbulence [5], nonlinear acoustics [6] 
and cosmology [7], etc. 

The nonlinearity of the KPZ equation includes the possibility of singularity formations in a finite 
time as a result of the local minima instability. Meaning that there is a competition between the diffusion 
smoothing effect ( the Laplacian term), and the enhancement of non-zero slopes. Let us mention the 
main properties of the KPZ equation for the cases v — > and v = 0. In the zero tension limit ( v — > 0) 
the KPZ equation has the following properties: (i) the unforced KPZ equation develops singularities for 
given dimensions. In one spatial dimension the sharp valleys are developed in a finite time t c . v ^o- The 
geometrical picture consists of a collection of sharp valleys intervening a series of hills in the stationary 
state [8]. In two spatial dimensions the KPZ equation develops three types of singularities in finite time. 
The first singularities are sharp valley lines with finite lengths, which the height gradients are discontinues 
while crossing the valley lines. The second type are the end points of the sharp valley lines. As time goes 
on these sharp valley lines hit each other and the crossing point of two valley lines produces a valley node. 
Generically these end points disappear at large time scales and only a network of sharp valley lines will 
survive [5,8,9]. In three and higher dimensions, the structure of the singularities can be more complicated. 
For instant, in three dimensions the singularities are, in the language of Burgers equation, shock surfaces, 
its boundaries, the intersection line of two shock surfaces and finally the point which three shock surfaces 
meet each other. The height gradients are discontinues while crossing the shock surfaces. A complete 
classification of the KPZ singularities, by considering the metamorphosis of singularities as time elapses, 
has been done in [10]. (ii) Similarly, for white in time and smooth in space forcing, in the zero tension 
limit singularities will be developed in a finite time in any given spatial dimensions. For instance, in 
two spatial dimensions the sharp valley lines are smooth curves where in the stationary state the sharp 
valley lines produces a curvilinear hexagonal lattice [9,15] (see Figs. (1)). In three dimensions it can be 
shown that the shape of the polyhedra tiling the space can not be determined uniquely, nevertheless, the 
minimal polyhedra in d = 3 will have 24 vertices [9] . 

The KPZ equation with vanishing surface tension ( v = ) produces multi-valued solution after 
time scale i C: „ = o [5]. In Fig. 2 we demonstrate the multi- valued solution of the unforced zero tension 
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KPZ equation in two dimensions. We have used the Lagrangian method to simulate the KPZ equation 
with v = and initial condition h(x,y,0) = sin(x)sin{y) [5]. The sinusoidal function is a typical initial 
condition and has been used only for simplicity. With this initial condition the time scale that the KPZ 
equation produces the multivalued solution can be found exactly as t c ,^=o = 1 [5]- As shown in Fig.(2-a), 
it is evident that for time scales t < t c „ =0 the height field is single- valued. At the time scale t — t c u=0 
the height field become singular (see Fig.(2-b)). Finally in Fig.(2-c) we have plotted the height field for 
time scale t > t c u=0 in which the height is multi-valued. The singularities in the limit v — > can be 
constructed from multi- valued solutions of the KPZ equation with v = by Maxwell cutting rule [5], 
which makes the discontinuity in the derivative of height field. Indeed the Maxwell cutting rule states, 
that for x's which the height field h(x, t) becomes multi-valued, the physical solution can be chosen so 
that, h p h(x, t) = Max{hi(x, t), /i2(x, t), ■ ■ ■ , /iat(x, i)}, where N is the number of the multiple solutions 
of the field h(x,t) at position x. 



L - 0.2 

1 = 0.1 




-I 



Figure 1: Different time snapshots of height fields of randomly driven two dimensional KPZ equation in 
the zero tension limit ( v — ► 0) . Two upper snapshots are belong to the height fields before singularity time 
scales t CiI/ _>o- In the lower figure which is for time scales t > t c ^ u ^o, the field h(x, y) is not diffcrentiable in 
singularity curves [12]. In this simulation the relation between the forcing length scale a and the sample 
size L, is a ~ L/3. The forcing strength D is equal to unity. 
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In this paper, an exact master equation is derived from the zero tension (y = 0) KPZ equation 
for the joint probability distribution function (PDF) of height and its gradients. The master equation 
enables us to determine the time evolution of the PDF of h — h and all of the moments < (h — h) n >. It 
is proved that the derived height's PDF for the v = case, shows lack of positivity after a finite time t c . 
Potential relation between the t c and the singularity time scale ( t c u ^ ) of the KPZ equation having an in- 
finitesimal surface tension is discussed. Details of calculations are presented in the appendices A,B and C. 



2 Joint PDF of height and its gradients for zero tension KPZ 
equation 

Let us define P(h, Ui,Pij, t) as the joint PDF of h = h — h, Ui = h Xi and Pij = h XiXj where i, j = 1, 2, • • • , d. 
Using the zero tension KPZ equation, it is shown in appendix A, that the P(h,Ui,Pij,t) satisfies the 
following equation 

Pt = l(t)P~ + f£ u 2 P~ a(d + 2)^ PllP 
I I 

-a ]T Plk p lm P Pkm + Ho)p~ - fc"(o) £ p UlUl 

l,k<m I 

+2fc"(0)£^ Pll +*r(0)£ P «*m 

/ Kk 



-2k""(0)J2P PllPkk - (4) 



Kk 

where "f(t) = h t , fc(x — x') = 2DqD(x. — x'), k (0) = k XiXi (0) and k (0) = k XiXiXiXi (0). Equation 
(4) enables us to calculate the joint PDF P(h, Ui,pij,t) and all the moments S n = (h n ) in rf-dimensions. 
Various moments has been calculated explicitly for the three-dimensional case in Appendix B. To derive 
a closed expression for P(h — h, Ui, t) one needs to know the moments such < h n u™-u l ^pij >. As appeared 
in Appendix B, using eq.(4) one can show that such moments are identically zero [11]. Using the identity 
mentioned above, it can be shown that P(h,Ui,t) has the following expression ( see appendix C); 



/d\ d d d 



(5) 



where, 



and, 



Z{\, m, t) = Fi(A, fii)F 2 (\, Hi)...F d {\, ii d ) exp {~X 2 k(0)t) (6) 



^•(A./ij.t) = (1 - t&nh 2 {^2ik xx {0)aXt))-i 



exp[- l -ak"(0)Xt 2 - ^i^ 2ik ^°\ & nh(^2ik xx (Q)aXt)]. 



(7) 



Indeed Z(X, Hi,t) is the generating function of the h and u^s. Therefore expanding eq.(6) in powers 
of A, all the moments S n can be derived. For instance, the first five moments are as follows 
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15 K ak"(0) n t 



(h 4 ) 




afc"(0) 




(A 5 ) 



x — d( d)(—) w + — d(— V] 

1 9 v 105 ,y t*' 3 t*' 1 



16 T/ 248 , w t Nin 32 „ i 



where t* = ( 



fc(0,0) 




The important content of the exact expressions derived above is that through them the time scale 
that the height's PDF P(h — h,t) lacks positivity condition after t c . The positivity of PDF means that 
all the even moments of < (h — h) 2n > must be positive. In fact the above moment relations indicate that 
different even order moments become negative in some distinct characteristic time scales. Closer looking 
at the even moment relations reveals that the higher the moments are, the smaller their characteristic 
time scales become such that asymptotically tends to t c = a^t* for very large even moments. The 
coefficients ad are of order of unity [8,11]. Indeed it can be shown that after time scale t c the right tail 
of the probability distribution function (PDF) of height fluctuations (i.e. P(h — h, t)) is going to become 
negative, which is reminiscent of the singularity creation. In what follows we argue that the two time 
scales t c and t c ^ =0 are related each other and t c ~ t CiV —o [8,11]. 

We note that the eq.(6) has the property that Z(0,0, 0,t) = 1 which means that P(h — 
h,u,v;t)d(h — h)dudv = 1 for every time t(0^t~<oo). So the PDF of h — h and its derivatives is 
always normalizable to unity. In the limit of v = after t CiV —o the height field becomes multi-value on 
the valleys, which is related to the left tail of the P(h — h). The multiplicity of height field on valleys, 
on which the height difference h — h is mostly negative, increases the probability measure in left tail 
of the PDF. Therefore to compensate the exceeded measure related to the multi-valued solutions the 
right tail of the PDF tails should become negative. Therefore one concludes that t c ~ t CtV -o. On the 
other hand as mentioned in the introduction the singularities in the limit v — > can be constructed from 
multi- valued solutions of the KPZ equation with v — by Maxwell cutting rule [5], which makes the 
discontinuity in the derivative of height field. Therefore the time scale that the zero tension KPZ equation 
produces multi-valued solutions is the same as the time scale of singularity formation in KPZ equation 
with infinitesimal surface tension. So t CyV =Q — t c ,i^o — t c , where t c — a^t* and t* = ( ^p^ ) ) 1 ^ 3 - The 
t c scales with a as a 3 . Hence the smaller the a, the shorter the time scale of singularity creation. 

Taking into account that a > and /c"(0,0) < 0, the odd order moments S^m+i are positive 
in time scales before formation of multivalued solution. It means that the probability density function 
P(h — h,t) in this time regime is positively skewed. Therefore the probability distribution functions of 
height difference has a non zero skewness as it evolves in time, at least up to the time scale where the 
multivalued solutions are formed. In Fig.(3-a), using the equations (4-7), we have numerically sketched 
the PDF evolution in time for zero tension 3 + 1-dimcnsional KPZ equation. As the system evolves in 
time, the formation of the multivalued solutions leads to the negativity of the right tail in the PDF. Also 
the evolution of PDF right tail for time scales before and after t c can be checked in Fig. (3-b) . 

Now we would like to add a few comments on the equation governing the PDF of h — h. As 
shown in eq.(4), for zero tension KPZ equation, the P(h,h x ,h xx ;t) satisfies a closed equation. Adding 
an infinitesimal surface tension ( v ^ 0) to the KPZ equation, the PDF equation will no longer be closed. 
Here we have proved that the height's PDF of zero tension KPZ equation lacks the positivity condition 
after the finite time scale t c . The situation is similar to the results given in [16], where the derived PDF 
of velocity increments of invisid Burgers equation is not positive definite. In 1+1 dimensions, it has been 
already proved that adding an infinitesimal surface tension will guarantee the positivity of the PDF of 



The same argument to find the singularity time scale i.e. t c0 , can be applied to the problem 



h-h [8,17]. 
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Figure 2: Different time snapshots of height fields of the zero tension {y = 0), unforced two dimensional 
KPZ equation with initial condition h(x,y,0) — sin(x)sin{y). In this case the time scale t c ,v=o is equal 
1. In the upper figure the height field is single valued ( for time scale t < t c , v =o)- In the middle we 
demonstrate the height field in time scale t = t c ,^=o which is the singularity time scale. In the lower 
figure, which is for time scales t > t c , v=0 , the field h(x,y) is multi-valued. 
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of decaying tensionless KPZ equation with random initial condition. We use the following probability 
density functional for initial height field (ho(x) and its spatial derivative (uq(x) in 1+1 dimension: 



PM*),».M] « =x P (-/WW. - *')M*') + /<W..(*W(* - *>.(*'/ 



(«) 



where the equality will be hold by a normalization constant and B"(x) = B xx {x). The initial distri- 
bution (8) shows that initial height field and its derivative are both Zero mean, statistically homogeneous, 
Gaussian processes and spatially correlated with covariance : 



<h (x)ha(x') >=2B(x-x') (9) 
<u (x)u (x')>=-2B"(x-x') (10) 
< h^(x)u^(x') >=< h%{x) >< u n (x') > . (11) 



We set B(x) to be Gaussian function with standard deviation <jq: 

(x 



B(x-x') = ^==eM- 2 )i (12) 



■Kok 



(To is the correlation length of initial height field. Now using the KPZ equation and its derivative 
we can derive the equations governing the time evolution of height and height derivative moments. This 
procedure has been represented in Appendix B ( by setting / = ) . Since we are interested in functional 
dependence of the time at which the singularities begin to form, it is sufficient to derive the second 
moment of height which can be easily found as: 

< h 2 > = -ia 2 < ul > 2 t 2 + < h a 2 > . (13) 

2 1/2 

Eq.(13) shows that the second height moment would be zero at time t$ = • Evaluat- 

ing < h > and < u 2 > from eq.(8) and using eq.(12) we see that scales as u\j 2 . Due to the fact that 
the time scale of singularity formation (i c o) is proportional to 4q by a constant of order of unity, the same 
conclusion is true for t C Q. Generalization of the above argument to higher dimensions is straightforward. 

d + 4 

It can be easily shown that t C Q ~ a 2 . This means that the smaller their characteristic initial length 
scales cause the singularities to be produced at smaller time scales. 




-4 -2 2 4 012345 
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Figure 3: In the upper graph the time evolution of PDF oih — h before formation of multivalued solutions 
at 0.2t c , 0At c and 0.8t c is numerically obtained. Lower graph shows the right tails of the PDF oih — h 
for 0.2£ c , 0At c and 1.2t c corresponding to time scales before and after formation of multivalued solutions. 
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To summarize, we obtain some results on the problem of KPZ equation in d + 1 dimensions with 
a Gaussian forcing, which is white in time and short-range-correlated in space. In the non-stationary 
regime, where the singularities are not yet developed, we find an exact form for the generating function of 
the joint fluctuations of height and height gradients. Starting from a flat initial condition, we determine 
the time scale of the singularity formation and the exact functional form of the time dependence in 
the height difference moments at any given order. We note that if there is no singularity in the height 
field, the time evolution of the height's PDF of tensionless KPZ and the KPZ equation with infinitesimal 
surface tension are identical. This is the reason that the term limit ^^(vV 2 h) is equal to zero for time 
scales before the singularity time scale t c . We were able to give the solution of the PDF's equation for 
time scales before t c . 

We believe that the analysis followed in this paper is also suitable for the zero temperature limit 
in the problem of directed polymer in a random potential with short range correlations [4] . 
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3 APPENDIX A 

In this appendix we prove the equation (4) . Define generating function Z(X, Hi, rjij, Xi, t) =< 0(A, fj,i,r]ij,Xi,t) > 
for the fields h = h — h, Ui = h Xi and pij = h XiXj . The A, in and rjij are the sources of h = h — h, ui = h Xi 
and = h XiXj , respectively and i, j = 1, 2, • • • , d. 
The explicit expression of 6 is as follows: 

d d 

6 = exp(-i\(h(x,y,z,t) — h(t)) - - i ^ VijPij) 

i— 1 i<j=l 

(14) 

where = h Xi and p^j = h XiXj . Also introduce qijk as = h XiXjXk which will be used later. 

By considering the zero-tension KPZ equation we will have the following time evolution for h and 
its derivatives; 

?H = fX>? + / (15) 

»=i 
d 

Ui, t = a^uipu + f Xi (16) 
i=i 

d d 

Pv-t = a ^2puPij + a^uiquj + f XiXj . (17) 
i=i i=i 

Also using the Novikov's theorem we can write the following identities; 

3 



(/6) = -*Afc(0, )Z-iJ2 mik"(0)Z (18) 
i=i 

(f Xi Q) = -iiHk"(0)Z (19) 

3 

</* 4 * 4 e> = -i\k"(p)z - 1 mik""(o)z (20) 

1=1 

(f XiXj e) = -ir ]ij k""(0)Z i^j (21) 

where 

k(x - x', y-y',z- z') = 2D D(x -x',y- y' , z - z') , 
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fc(0) = fc(0,0,0) = -^ 
k xxxx (0,0,0) = f^o. 



fc'(O) = fc, (0, 0, 0) = k y (0, 0, 0) = k z {0, 0, 0) = 0. 

Differentiate the generating function Z with respect to t and using the eqs.(9-f 5) and the following 
identity 

iZ Xl - i\Z H - i 53 ViZmi = 51 7 fe (<Hfl®)i ( 22 ) 

i i<j 

time evolution of Z can be written simply as 

Z t = iXj{t)Z - i^-^Z^ -ia^Z m +ia^ V*jZ mm - X 2 k(0)Z + 53 tfk"(0)Z 

l l Li<j I 

- 2\J2mik"(0)z - C£mivkk + J2rfk)k""(0)z, (23) 

l l,k Kk 

where 7(f) is defined as 7(f) = h t . Fourier transforming Z respect to X,fii and 77^, the joint probability 
density function (PDF) of h,Ui,pij, P(h,Ui,Pij,t) is generated as; 

P(h,Ui,pij,t) = [ ^ TT ^ TT exp(iA(/i(a;, y, z, t) +h(t)) +i^fJ>iUi+iY2 r )ikPik) 

/ Z7T Z7T Z7T ' — ' ' — ' 

17 t i<j I l<k 

xZ(X,fii,i]ij,Xi,t) (24) 

By the use of equations (|23|) and 124|l the equation governing the P(h,Ui,pij,tys time evolution will be 
derived as 

Pt = l(t)P~ + I ^ _ a(d + 2 ) £ PlI P -a^ PikPi m P Pkm + k(0)P~ - fc"(0) J] 

i I l,k<.rn I 

+2k"(0) ^ P~ pu + fc""(0) ^ P PifcPifc - 2fc""(0) 53 P piIJJM (25) 

I l<k Kk 

where d is the spatial dimension. 

Eq. (|25|l enables one to write the time evolution of the moments of height and it's corresponding 
derivatives. The obtained equation has the following form 

f (h n °AB) = -THrrWGfr^AB) - ^53<ft"°- 1 J 4I?Kf>+a 53 n lm (/i"M/*) 

I l,k<m Pkm 

,h n °AB> 



-a 
I 



^2{h n °AB Pll )+k(0)n (n a - \){h n °- 2 AB) - fc"(0)5>i(n, - 1)( — } 



+2k"(0) 53 n n„( — ) + fc""(0) 53 n lk (n lk - 1)( — ) 

1 p « ;<fe 

^ f.n a A r> 

+2k""(0) 53 n u n kk ( ) (26) 

ttk mPkk 



where 



A = U U " 

i=l 

By choosing different sort of no, 7ij and mj values, various type of coupled differential equations, 
governing the evolution of the moments, can be constructed . 
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4 Appendix B 



In appendix A we have derived the general time evolution equation for the moments for the arbitrary 
(d+l)-dimensional case. The moments in 1+1 and 2+1 dimensions has been derived in refs. [8] and [11], 
respectively. Here we are going to restrict ourselves to the (3+l)-dimensional case. Also the moments in 
d+1 dimensions is given at the end of this appendix. 

In this appendix the height moments will be calculated exactly and it will be shown that (h XiXj 0} , i + 1 
j will be zero by considering flat initial condition in the (3+l)-dimensional case. However it can be shown 
that in any general (d+l)-dimensional case this identity will also be true. By choosing different sort of no, 
rii and n^- values, various type of coupled differential equations, governing the evolution of the moments, 
can be constructed. The first example is to choose hq = m = Uij = then; 

-aJ2(Pu) = -"(V-u) =0^ (v-u) = (27) 
i 

which verifies that the fluid is incompressible in average. The main aim is to calculate the moments 
(Pij exp(—i\h — i J2i Mz u z))- Before that, we have to follow several steps. First of all we have to calculate 
the moments such as {pijUiUj) and i ^ j. Using the eq.(20), we have 

4; a, a , P,. = a^iuiUjpupij) - a^2{uiUjPuPij) = a'^(u i u.j(pupi j -puPij).) (28) 
l i i 

Looking at the right hand of equation (|28|l we see that the terms with / = i or I = j cancel each 
other. Noting that we have restricted ourselves to the 3+1 dimension, the equation (|28|l can be written 

^(uiUjPij) = (u.Ujipiipij - puPij)) i^j=£l (29) 
Now inserting the rhs of equation l|29|) in equation (|26|l . one finds 

d 
di 

a^iuiUjPkiPkipij) + a ^(UjUjPkjPkipn) 
k k 

-ay^UjUjPkkPuPij) - a/^UjUjPkiPkiPij) 

k k 

-a^juiUjPkiPkjPu) + a ^(ujUjPkkPtjPu) (30) 

k k 

It can be easily seen that in the case that i + 1 j + 1 I the rhs of equation (j30|l will be zero. Using 
this result and considering a flat initial condition we have; 

{uiUjipupij -puPij)) = (31) 

and 

(iitUjPtj) = 0. (32) 

Also it can be shown by induction that all the moments such as (h no p™j 3 u^u™ 3 ) are zero too. 

Therefore one concludes that (p^ exp (—i\h — i ^ p-iui)} = 0. As shown in appendix C this relation is 
crucial to derive the eq.(6). 

Now let us calculate the moments of u^'s. Using the eq.(2) it can be shown that 

!«) = -a(u?X>«> - n(n l)fc"(0)«- 2 ). (33) 



-(UiUjipiiPlj - PuPij)} 
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Differentiating and (ufuj) with respect to Xi and Xj and using the statistical homogeneity 

and equation (26) we have 

(u?p«) = (u?p«) = 0. (34) 

Therepore; 

|(^) = -n i (n i -l)fc"(0)« i - 2 ) (35) 

which is an iterative equation implying that any order of the Ui moment can be calculated by 
knowing the lower moments. Because (ui) — 0, from equation (|35|l it is obvious that any odd moment of 
Ui will be zero. For rn — 2,4,6,8, (u™*) will be 

(u 2 ) = -2fc"(0)t (36) 

{ui) = 12fc"(0) 2 t 2 (37) 

(«?) = -120fc"(0) 3 t 3 (38) 

(«?) = 1680fc"(0) 4 t 4 . (39) 

and 7(i) ( 7 (t) = 7i(f)) will be; 

7(t) = ^{uf) = -3ak"(0)t (40) 
For moments such as (u^u™ 3 ), it will be deduced that 

§( u >?) = -«i(«< - i)fc"(o)« <_ V> 

—2 



- nj (n,- - l)fc"(0)«'^- 2 ) - a«*t^ $>»> (41) 

i 

Differentiate (u^u™ 3 ) and (u^u™ 3 Uk) with respect to £j and a;*,, respectively one finds 

= n i « i - 1 ^p ii )+ r^^-Vj} W (42) 



{u i % u j :l Uk) Xk = ni{u i * Uj'ukPik) +n j (u i t u :j 3 UkPjk) 

+ {uTu n J 3 Pkk ) i^j^k. (43) 

Using the statistical homogeneity and identity such as {uiUjPu) = ( I ^ i,j) it can be seen that 
equation i|4*T|) would be 

= - l)fc"(0)<C"V> 

-n^nj-lJA/'COXu^u^- 2 ) (44) 

Therefore one finds the following expression for the moments; 



(u 2 u 2 ) = 4k" 2 (0)t 2 (45) 
(u 4 m 2 ) = -24fc" 3 (0)t 3 (46) 
(ufu 2 ) = 240fc" 4 (0)t 4 (47) 



and 
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(nfuiui) = -8k"\0)t\ (48) 

Now we want to calculate (ft 2 ) and its higher moments. Using eq.(20) it can be shown that the 
moment (ft 2 ) satisfies the following equation 

j t (ft 2 ) = -a^Chuf) -aJ2(h 2 pil) + 2fe(0) (49) 

But the moment (h 2 pu} can be written as 

(ftV) = (h 2 Ul ) Xl - 2{huf) (50) 

where (h 2 ui) Xl = 0, so 

(h 2 Pll ) = -2(hu 2 ) (51) 
Then equation (|49|l can be written as 

j t (h 2 )=aJ2(hu?) + 2k(0) (52) 
Before calculating the (huf) moment, the more general term (h^u™ 1 ), can be studied, so 

j t (h na uT) = -n^m^u?) - ^^(ft- M - u 2 ) 

-aJ2{h n °uTPii)+n (n - l)k(0)(h n °~ 2 uT) 
i 

-mim - i)fc"(o)(ft"°<'- 2 ) (53) 

where 

Ch n °uTPii) = -^-(h n °- 1 ur +2 ) (54) 
(h no u? Pjj ) = -n^-'uTu]) i+j (55) 
Therefore equation (|53[) can be written in a new form as 

j t (h n °u?) = -n o7 (t)(ft"- 1 <') + ^ ^{h^uTu 2 ) 



ano 



^^(ft' io - 1 7i^ U 2 )+n (n_l)fc(0)(ft"«- 2 u^) 



-n t { ni - l)fc"(0)<ft"»<*- 2 ) (56) 
From equation lfS^|l . {huf) can be easily obtained. 

4 £(ftu 2 ) = -y(t)J>?> 



i 



6 

So finally by using equations dlt and l|45(l one gets; 
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J2( hu h = ~4ak" 2 (0)t 3 (58) 
l 

Therefore we find; 

(h 2 ) = -a 2 k" 2 (0)t 4 + 2k(0)t (59) 
Now it can be shown that the moment (h 3 ) satisfies the following equation; 

j(h 3 ) = -3j(t)(h 2 ) 3 | Y,Ch 2 uf) a^(h 3 Pu) (60) 

l i 

By spatial differentiation (h 3 pu) will be proportional to (h 2 uf), so 

|(P) = -3 7 (t)(P)+3|^(Pn?) (61) 

i 

and in a similar way the time evolution of moment (h 4 ) can be written as 
d 



dt 



(h 4 ) = -^(t)(h 3 ) + 2aJ2(h 2 u 2 ) + 12k(0)(h 2 ). (62) 



l 



It is appear that for calculating the moments (h 3 ) and (h 4 ), we should calculate the moments 
(h 3 uf) and {h 2 uf}. For (h 3 ) it is necessary to know ^i{h 2 u 2 ), so 



|E^V> = -27(i)£</^> - f £<*«f> 
+2a£(/n« + 2fc(0)^>f) - 6fc"(0)£(P) (63) 



Kk I I 

where (huf) and (hufu 2 ) should be calculated from the following equations 



|E(^) = -7(*)E^)-^E(^) 



;<fc ;<fc ;<fe 

+^<«?«i«l>-4fc"(o)E^«?>- 



(65) 



Using the eq.(52) we find; 



(/i 3 > = ~a 3 k" 3 (0)t e . (66) 
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To calculate the moment (h 4 ) one needs the following moments 



(feui«2«3> = -16a/c" 4 (0)i 5 (67) 
^2(huf) = -720ak"\0)t 5 (68) 



5>M)=432fc"¥ (69) 

5^<«f«|w|u?) = 144fc" 4 t 4 (70) 
i 

£<« = 1440fc" 4 t 4 (71) 

^ (huful) = -288k" 4 t 5 (72) 



J2(h 2 uf) = ^ a 2 k" 4 {0)t 6 + 72/c(0)fc" 2 (0)t 3 (73) 
i 

y2(h 2 ufu 2 k ) = ^-a 2 k" 4 t 6 + 48k(0)k" 2 (0)t 3 (74) 
z<fc 

J2(h 3 u 2 ) = ^a 3 k"H r - 24afc(0)fc" 2 (0)t 4 (75) 
so, finally by substituting the moments above in eauation (|62|l we find; 

(h 4 ) = ~—a 4 k" 4 t s - 12a 2 k(0)k" 2 (0)t 5 + 12k(0) 2 t 2 (76) 
35 

Generalizing this method to the d + 1— dimensional case by a similar amount of calculations, the 
second, third and fourth moments can be written as; 

where t, = ( Jfii^;)'". 

5 Appendix C 

In this appendix using the identities which have found in appendix B, joint-probability distribution 
function (PDF) is calculated for zero tension KPZ equation. Similar to appendix B we restrict ourselves 
to the 3+1 dimensions case. 

The zero tension KPZ equation in 3+1 dimensions has the following form; 

ht ( x ,y,z,t)-~(hl + hl, + h 2 z ) = f (77) 

Now defining 
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h x = u,h y = v, h z = w (78) 
Differentiating the KPZ equation (|77() with respect to x,y and z, we have 

h t = %(hl + h 2 v + hl)+f(x,y,t) (79) 
u t = a(uu x + vv x + ww x ) + f x (80) 

Vf = a(ltUy + VVy + WWy) + fy (8l) 

w t = a(uu z + vv z +ww z ) + f z (82) 

for h and corresponding velocity fields. The generating function Z(X, fj,i, fi%, x, y, z, t) is defined 
such as to generate the height and velocity field moments. By introducing as 

9 = exp(—i\(h(x, y, z, t) — h(t)) — i[i\u — i^v — ifXaw) (83) 

The generating function will be written as Z(X, /ii, A*3i x i z i t) = (©)■ 

Using the KPZ equation (|77?f and its differentiations (78), (79) and (80) with respect to x , y and 
z, the time evolution of Z(A, fii, fJ-2, M3j x i Hi z i can be written as 

Z t = i"f{t)\Z — i\—((u 2 + v 2 + w 2 )Q) — ia(ix((uu x + vv x + ww x )Q) — ia^Hwu-y + vv y + ww y )Q) 

-ia^((uu z + vv z + ww z )Q) - iX(f&) - ifii(f x &) - i^ify®) - i^ifz®) (84) 

where "f(t) — h t = j(u 2 + v 2 + w 2 ). By considering statistical homogeneity we have 

Z x — ({—i\u — i[i\u x — i[i 2 v x — i[i3w x )@) (85) 
Z y = ((—i\v — ifiiUy — i[i 2 v y — in$Wy)&) (86) 
Z z = {(-iXw - i\i\u z - i/-J. 2 v z - ifi3W z )Q) (87) 

Because we are in a time regime that the singularities has not been formed yet, the order of partial 
derivatives can be exchanged 

Qx-dx = dx dx so Vx ~ u y> w x ~ u z an d w y — v z- Keeping the definition of l|83|) in mind, we 
can easily write 

d 

i- — {(-iniu x - i(j, 2 v x - ifi 3 w x )Q) = (u x Q) - 
dm 

i/j,i(uu x Q) - i/x 2 (ww K 0) - i(j, 3 (uw x ®}. (88) 
From equations (|85|l and Ij88(l we have 

(u x &) - i^i(uu x Q) - ifi 2 (uv x <3) - i^, 3 (uw x &) 

= -\—(uO) = -i\Z^ 1 (89) 

And in a similar manner one finds; 

- iAZ M2M2 = (v y Q) - ifxi(vu y @) - ifi 2 (vv y Q) - ip, 3 (vw x Q) (90) 

- iAZ M3/l3 = (w z Q) - ifii{wu z O) - ipL 2 (wv z Q) - i^i 3 (ww z O) (91) 

By using the Novikov's theorem the expression (/©} and (f Xi &) can be written respect to Z ( 
appeared in appendix A ). So we have 

a d ot 

Z t = ij(t)\Z — i\—Z flltll — i^—Z^ 2ll2 — iX—Z^^ — a(u x Q) — a(v y <d) — a(w z Q) 

- \ 2 k(0,0,0)Z + t i 2 1 k xx (0,0,0)Z + ^ 2 k xx {0,0,0)Z + n 2 3 k xx {0,0,0)Z. (92) 
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The (u Xi Q) terms can be written as 



(u x Q) = — {Q) x H ((iXu + i/j,2V x + i[i 3 w x )Q) 

Mi Mi 

- -i—Z^-^(v x e)-^(w x &), (93) 
Mi Mi Mi 



and similarly for (v y @) and (w z @) we have 

X ry Mi / ~\ M3 



(tv8) = -i-^-q«,e>-fKe) (94) 

M2 M2 M2 



(^e) = _iA^ ^_a< u ,e) - ^-(v z e) (95) 

M3 M3 M3 

The terms such as (h XiXj Q) (i ^ j) are the main troublesome terms which appear in the Z l s time 
evolution equation ( i.e (u y Q) ), preventing us to write the Z— equation in a closed form. Fortunately as 
it is shown in appendix B, these terms will become zero by considering a flat initial condition. 

Therefore Z satisfies the following equation; 

Zt = ij(t)XZ - iA^Z MlMl - iX^-Z^^ - iX^Z^ 3fl3 + ia—Z^ - ia — Z^ 2 
i i z Mi M2 

~ia—Z M ~ A 2 fc(0, 0, 0)Z + (4k xx (0, 0, 0)Z + f4k xx (0, 0, 0)Z + (4k xx (0, 0, 0)Z (96) 

M3 

In what follows we are going to solve the partial differential equation above, by using a flat initial 
condition, h(x, y, z, 0) = u(x, y, z, 0) = v(x, y, z, 0) = w{x, y, z, 0) = 0, which equivalently means to write 

P(h, u, v, 0) = 5(h)6{u)5(v)6(w). (97) 

This means that; 



2(0,0,0,*) = 1 (98) 

A useful and efficient way to solve the Z, time-evolution differential equation is to factorize it in 
the following manner [8,11]; 

Z(X,(j,i,H2,H3,t) = Fi(X, fi 1 ,t)F 2 (X, n2,t)F 3 (X, /z 3 ,*) 

xexp(-A 2 fc(0)i) (99) 

Then by inserting (|99|) in (|96|l we obtain 



■Fit -Fiji's + FiFitFs + F\F2F 3t — 

^F 2 F 1 ^ 1 ^ 1 F 2 F 3 - iX- 
a „ „ „ A „ „ „ A 



i 1 {t)XF l F 2 F 3 - iX-F 2 F llMllJll F2F 3 - iX-F^^Fs 



-iX-F^Fs + ia—F lhll F 2 F 3 - ia—F 1 F 2lJ , 2 F 3 
l Mi M2 



Aa—FxFiF^ - X 2 k{0)F 1 F 2 F 2 + l i 2 1 k"(0,0)F 1 F 2 F 3 



X_ 

M3 

+fi 2 k"{0)F 1 F 2 F 2 + ^k"{Q)F 1 F 2 F 2 . (100) 



Therefore; 
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F t = -iX-F^ + ia-Fp + [fJ?k"(0) - iaXk"{0)t]F 

(101) 

with the initial condition F(A,/x,0) = 1. This points out that the height gradients in the three 
dimensions evolve separately from each other before the shock formations, and they are only coupled 
with the height field. By Fourier transforming equation l|l(-)l[) respect to [i a simpler partial differential 
equation of order one will appear. This equation can be solved by the method of Characteristics [8]. 
Finally the solution of F will be 

F(fi,X,t) = (1 - tanh 2 (y/2ik xx (0)aXt))~i 

ex P [-^"(0)Ai 2 - ^J^M tanH ^/ 2 ik xx (0)aXt)} 
2 2 V u\ 

(102) 

and 



Z(X,/j,i,H2,t) 

= F(X, m,t)F(\, » 2 ,t)F(X, fi 3 ,t) exp (-A 2 fc(0, 0)t) 
= (1 - t&n\i 2 {^2ik xx {Q)aXt)y i 



exp[~ j-^ 2 + /i 2 + f il)J 2lkxx } h aiih{^/2tk xx (0)aXt)} 
2 V a\ 

exp[-yafc"(0)Ai 2 - fc(0)A 2 t] (103) 

By inverse Fourier transformation of the generating function Z, the probability distribution func- 
tion (PDF) of the height fluctuation can be easily derived; 

vru 4-\ f dX dni dfi 2 dn 3 7 . 

r{n,u,v,w,t) = / exp(zA/i + i[i\u + t^v + i^w)Z 

J 2tt 2tt 2it 2h 

(104) 

Expanding the solution of the generating function in powers of A, all the {(h — h) n ) moments can 
be derived. For instance, the first five moments before the sharp valley formations are 

<S 2 >=^>)(J^[-(±)< + 2±] (105) 
x[-|(^-) 8 -12(^) 5 + 12(p) 2 ] (107) 

« = -« 5 *" 5 «»( gggfcj)* 

*hnF<F> 10+32 <F> 7 i < 108 > 



where t - ( fc(0 ' 0) 'l 1 ^ 3 
wneie - (, Q 2 fc /, 2(0 ^ ) 
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